SPATIAL DISTRIBUTION OF VOLCANIC HOTSPOTS AND PATERAE ON IO: IMPLICATIONS FOR 
TIDAL HEATING MODELS AND MAGMATIC PATHWAYS. C. W. Hamilton 1 , C. D. Beggan 2 , R. Lopes 3 , 
D. A. Williams 4 , and J. Radebaugh 5 , 'Planetary Geodynamics Laboratory, NASA Goddard Space Flight Center, 
Greenbelt, MD, 20771, USA, hamilton@daI.ca, 2 British Geological Survey, Edinburgh, UK, 3 Jet Propulsion Labor- 
atory/Caltech, Pasadena, USA, 4 School of Earth and Space Exploration, Arizona State University, Tempe, USA, 
5 Geological Sciences, Brigham Young University, Provo, USA. 


Introduction: Io, the innermost of Jupiter's Gali- 
lean satellites, is the most volcanically active body in 
the Solar System. Io's global mean heat flow is ~2 W 
m' 2 [1], which is ~20 times larger than on Earth. High 
surface temperatures concentrate within "hotspots" 
and, to date, 172 Ionian hotspots have been identified 
by spacecraft and Earth-based telescopes [2]. The La- 
place resonance between Io, Europa, and Ganymede 
maintains these satellites in noncircular orbits and 
causes displacement of their tidal bulges as the over- 
head position of Jupiter changes for each moon [3, 4], 
Gravitational interactions between Jupiter and Io dom- 
inate the orbital evolution of the Laplacian system and 
generate enormous heat within Io as tidal energy is 
dissipated. If this energy were transferred out of Io at 
the same rate as it is generated, then the associated 
surface heat flux would be 2.24 ±0.45W nf 2 [5]. This 
estimate is in good agreement with observed global 
heat flow [e.g., 1 ], but to better constrain tidal dissipa- 
tion mechanisms and infer how thermal energy is 
transferred to Io’s surface, it is critical to closely ex- 
amine the spatial distribution of volcanic features. 

End-member tidal dissipation models either consid- 
er that heating occurs completely in the mantle, or 
completely in the asthenosphere [6]. Mixed models 
typically favor one-third mantle and two-thirds asthen- 
osphere heating [7], Recent models also consider the 
effects of mant 1 e-as then o s ph ere boundary permeabil- 
ity and asthenospheric instabilities [4], Deep-mantle 
heating models predict maximum surface heat flux 
near the poles, whereas asthenosphere heating models 
predict maxima near the equator — particularly in the 
Sub-Jovian and Anti-Jovian hemispheres, with smaller 
maxima occurring at orbit tangent longitudes [4], Pre- 
vious studies [e.g., 1, 4, 8] have examined the global 
distribution of Ionian hotspots and patera (i.e., irregu- 
lar or complex craters with scalloped edges that are 
generally interpreted to be volcanic calderas [9]), but 
in this study, we combine a new geospatial analysis 
technique with an improved hotspot and paterae data- 
base [2], 

Methodology: Our geospatial analysis calculates 
nearest neighbor (NN) distances along great circles 
projected onto the surface of spheres or ellipsoids that 
correspond to the geometry of each of the major plane- 
tary bodies within the Solar System. In this case, we 
use an ellipsoidal geometry for Io, with a semi-major 
axis of 1829.7 km and semi-minor axis of 1819.2 km. 


Our program estimates the expected mean NN dis- 
tance for a given sample-size and region. The ratio (R) 
between the measured mean NN distance (r a ) and ex- 
pected mean NN distance (r e ) is used in combination 
with the standard deviations of r e to assess if the inputs 
are consistent with the an expected distribution model 
(r a - r e ), clustered with respect to the model (r a « ;* e ), 
or repelled relative to the model (r a » r e ). A second 
test statistic, c, evaluates the significance of the result 
implied by R. To calculate c, the difference between r a 
and r e is divided by the expected standard error [10, 
11]. A variety of expected models of spatial organiza- 
tion can be tested, but we focus on Poisson (i.e., ran- 
dom) distributions. For a given population size and 
region, if c is in its ±2c limits and R is in its ±2g lim- 
its, then the input distribution fits the expected model, 
whereas if c and R are outside their respective ±2a 
limits, then the input distribution exhibits a significant 
departure from the model. In such cases, R can either 
be less than -2a, implying smaller than expected (i.e., 
clustered) NN distances, or R can be greater than +2a, 
implying that the NN distances are greater than ex- 
pected (i.e., repelled). If c is in its ±2g limits and R 
outside its ±2 g limits, or c is outside its ±2a limits and 
R in its ±2o limits, then the results are not significant 
at the 2a level. We apply this NN technique to 172 
hotspots and 529 paterae on Io [2] to examine global 
and hemispheric distributions of volcanic features in 
each of the following domains: (1) Global, (2) Sub- 
Jovian hemisphere, (3) Anti-Jovian hemisphere, (4) 
Leading hemisphere, (5) Trailing hemisphere, (6) 
Northern hemisphere, and (7) Southern hemisphere. 
This allows us to test tidal dissipation models by de- 
termining if volcanic features on Io are randomly dis- 
tributed, or if they exhibit predicted patterns of spatial 
organization that are related to magmatic upwelling. 

Results: Hotspots and paterae have modal NN 
distances of 200-300 km and <50 km, respectively. 
Results for R (Fig. 1) and c (Fig. 2) show that the 
global distribution of hotspots, and hotspots in each of 
the six hemispheres, have NN distributions that are 
consistent with Poisson (i.e., random) models within 
la thresholds of significance in R and c (n.b., a fit 
within ±la implies a stronger fit than one within ±2o 
thresholds). In contrast, paterae have NN distances that 
are clustered relative to a Poisson model. 
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Figure 1. R-values of paterae (P.) and hotspots (H.) on lo. In 
Figs. 1 and 2, the black lines show ideal values for R and c, 
respectively, given a range of sample-sizes, whereas the dark 
grey and light grey lines show the ±1 a and ±2a thresholds. 
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Figure 2. c-values of paterae (P.) and hotspots (H.) on lo. To 
fit the model of a Poisson NN distribution both R- and c-values 
must be within thier respective ±2a limits. 


Unlike hotspots, paterae exhibit non-random distri- 
butions that are significant beyond 2o thresholds for R 
and c in all regions, except for the Sub-Jovian hemi- 
sphere. In the Sub-Jovian hemisphere, paterae have 
NN distances that are described by a Poisson model 
within 1 and 2o thresholds of significance in R and c. 

Discussion: Observations from Voyager and Gali- 
leo have been used to show that Ionian hotspots appear 
to have a "uniform" distribution that is not consistent 
with mantle dominated heating models [1, 8], Howev- 
er, we caution that these previous studies used "uni- 
form" to describe a uniform random process, rather 
than a uniform distribution in the sense used by [11], 
According to [11], a uniform distribution would imply 
strong repelling because all points would exhibit equal 
NN separations. This is not what was meant by [1, 8], 
who instead used the word "uniform" to mean "Pois- 
son". Our results show that globally, and within all 
hemispheric domains, hotspot locations are consistent 
with Poisson distribution models. This agrees with [1, 


8] and supports tidal dissipation models involving as- 
thenosphere or asthenosphere-dominated heating [e.g., 
4, 12], but not mantle-dominated heating models that 
would give rise to polar clustering of hotspots. 

Interpretation: Hotspots provide evidence of cur- 
rent volcanic activity, whereas paterae include volcan- 
ic surface features formed over longer periods of geo- 
logic history. Given that hotspot locations are con- 
sistent with Poisson distribution models, and paterae 
show pronounced clustering, hotspots may be fed by 
randomly distributed magmatic upwelling zones, typi- 
cally separated by 200-300 km, whereas paterae clus- 
ters could form over longer timescales as upwelling 
magma pathways branch out over lateral distances <50 
km. 

Spherical harmonic analysis of the global distribu- 
tion of mountains and volcanic centers on lo [13] 
shows statistically significant power at degree 2 with 
anti-correlations between mountains and volcanic cen- 
ters. The correlations between mountains and volcanic 
center distributions at high spectral degree could be 
related to structural links between some mountain 
blocks and magmatic systems, whereas the anti- 
correlation at low degree implies that most volcanic 
features form independently of the mountains. The 
independent formation of volcanic centers is consistent 
with our observation that mean NN distances between 
hotspots appear random. However, structural controls 
in the vicinity of magma upwelling zones could ex- 
plain clustering of paterae into preferred volcanic sub- 
domains over longer timescales. Global geologic map- 
ping supports the theory that hotspots concentrate 
within patera-forming regions [14]. 

Conclusions: Nearest neighbor (NN) analysis 
shows that Ionian hotspot distributions are consistent 
with Poisson (i.e., random) models, whereas paterae 
generally appear clustered relative to Poisson distribu- 
tions. These results support asthenospheric, or domi- 
nantly asthenospheric heating models, and suggest that 
while locations of active volcanism (i.e., hotspots) 
appear randomly distributed, patera regions may de- 
velop clustering as upwelling magma pathways branch 
out at shallow depths to feed new volcanic centers. 
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